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Abstract 


This  work  reviews  the  literature  on  gaseous  flow  pressure  resis¬ 
tance  in  packed  beds  and  found  important  differences  depending  on  test 
conditions  and  Reynolds  number  ranges.  A  test  apparatus  was  constructed 
which  allowed  for  the  testing  over  a  wide  range  of  pressures,  test  condi¬ 
tions,  and  Reynolds  number  range  several  orders  of  magnitude  higher,  than 
previously  tested.  From  the  resulting  data  it  was  ascertained  that  the 
classical  Reynolds  number  dependency  of  the  coefficient  of  drag  is  not 
correct  for  Reynolds  numbers  greater  than  10^.  A  new  correlation  for 
the  coefficient  of  drag,  given  below,  was  developed,  namely, 

Fv  =  3.33  x  106'DU  Re/U2 

b  avg 
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This  formula  was  shown  to  be  valid  for  a  Reynolds  number  range  10  -10  , 
for  spherical  particles  ranging  in  diameter  from  0.96mm  through  6mm. 

This  correlation  is  expected  to  provide  better  input  for  more  accurate 
calculations  of  transient  two-phase  flows  at  high  pressures. 


Mention  o{  the  commercial  products  used  in  connection  ukth  the  toork 
reported  here  do  cm  not  constitute  an  endors  ement  by  the  University 
o &  Illinois  or  the  contracting  agency. 
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CHAPTER  1 

INTRODUCTION 


The  gaseous  flow  resistance  through  packed  beds  of  various  sized 
solid  particles  has  been  investigated  by  numerous  researchers,  both  theo¬ 
retically  and  experimentally.  This  work  extends  back  to  the  days  of 
Osborne  Reynolds  in  the  early  1900's.  The  factors  which  determine  the 
energy  loss  (or  pressure  drop)  are  numerous,  and  the  analysis  require 
many  simplifying  assumptions.  However,  the  four  main  factors  that  directly 
affect  the  resistance  of  gaseous  flow  through  a  packed  bed  are: 

1.  Rate  of  fluid  flow 

2.  Viscosity  and  density  of  the  fluid 

3.  Porosity  and  orientation  of  packing 

4.  Size,  shape  and  surface  roughness  of  particles 

The  first  two  factors  pertain  to  the  fluid  and  the  last  two  the  solids. 
Fluid  Properties 

The  pressure  drop  through  a  granular  bed  is  proportional  to  the  fluid 
velocity  at  low  flow  rates  and  approximately  to  the  square  of  the  velocity 
at  higher  flow  rates.  Reynolds  was  the  first  to  derive  a  formula  that  re¬ 
lated  fluid  resistance  (pressure  gradient)  to  motion,  due  to  friction  in 
the  form: 

AP/L  =  au  +  bpu2  (1) 

where  AP  is  a  pressure  drop  over  a  length  L,  p  is  the  density  of  the  fluid, 
u  an  average  gas  velocity,  and  £  and  b.  are  constants.  Dividing  this  equa¬ 
tion  by  velocity  yields: 

AP/Lu  =  a  +  bpu  (la) 
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A  plot  of  (AP/Lu)  versus  (pu)  produces  a  straight  line.  It  was  found 
that  this  expression  accurately  represented  the  relation  between  flow  rate 
and  pressure  drop  at  very  low  Reynolds  numbers.  As  the  velocity  approaches 
zero,  it  is  seen  that  the  ratio  of  pressure  drop  to  velocity  is  a  constant. 

lim  AP/L 

.  — - —  =  a 

u  -*■  0  u 


According  to  Darcy's  law,  this  contant  is  proportional  to  viscosity  (p) . 
Taking  the  limit  at  high  flow  rates,  the  constant  a  is  negligible  com¬ 
pared  to  bpu,  where  pu  is  the  mass  flow  rate. 


lim 


AP/L 


=  bpu 


(lb) 


This  condition  exists  in  turbulent  flow  where  the  flow  resistance  is  due 
to  kinetic  energy  losses.  Hence,  the  above  equation  can  be  rewritten  as 
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AP/L  =  a  pu  +  bpu 


(lc) 


where  a'  -  a/p  (a  factor  pertaining  to  variables  in  solids  only).  The 
first  term  on  the  right  hand  side  represents  viscous  energy  losses  and 
the  last  term,  kinetic  energy  losses. 


Solid  Properties 

Blake  was  the  first  (1914)  to  establish  the  dependency  of  pressure 
drop  upon  porosity,  by  using  an  analogy  to  a  circular  pipe.  He  developed 

two  dimensionless  groups: 

AP  °h  j.3  D  pu 

Group  A  =  — j  *  i  Group  B  =  (Id) 

pu  J 

Modified  Friction  Factor  Modified  Reynolds  Number 


Here  porosity,  <J>,  represents  volume  fraction  for  the  gases  to  penetrate 


3 


into  the  packed  bed.  See  Figure  la,  for  schematic  of  a  packed  bed  and 
explanation  of  terms.  Note  that. 


Volume  Solids 
Volume  Chamber 


(2) 


where  N  =  number  of  beads  in  the  volume  bounded  by  the  cross  sectional 
2 

area  (ttDc  /4)  and  length  L. 

The  calculated  porosity  is  only  an  average  porosity  for  the  packed 
bed.  Figure  lb,  from  Benenati  and  Brosilow,  shows  how  the  void  fraction 
varies  through  a  packed  bed  [2].  Since  the  sphere  only  makes  point  con¬ 
tact  with  the  wall,  the  porosity  approaches  1  at  the  wall  and  then  oscil¬ 
lates  until  it  dampens  out  at  the  average  porosity  of  the  packed  bed,  about 
4  1/2  to  5  sphere  diameters  distance  away  from  the  test  section  wall. 

Blake  obtained  a  relation  between  Group  A  and  Group  B  by  plotting 
one  against  the  other.  His  efforts  failed  because  he  did  not  realize 
the  pressure  drop  was  due  to  both  kinetic  and  viscous  energy  losses 
simultaneously  [1]. 

Later  Burke  and  Plummer  theorized  that  the  total  resistance  was  equal 
to  the  sum  of  individual  particle  resistance.  They  stated  that  viscous 
energy  loss  was  proportional  to  ^  and  kinetic  energy  loss  to  . 

<T  <f> 

However,  they  failed  to  realize  the  additive  affects  of  these  losses  [1]. 

2 

Kozeny  arrived  at  viscous  energy  loss  proportional  to  ,  by 

<T 

assuming  that  the  granular  bed  was  equivalent  to  a  group  of  similar  paral¬ 
lel  channels. 

In  1949,  Ergun  and  Orning  furthered  the  work  in  this  area  and  derived 


Void  Fraction  » 
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>n  equation  for  drag  per  unit  volume  as  follows  [3]: 


A]1  -  k-  .  i-LiT  . 

L  1  ,3 

4> 


+  K  •  - I-Z-  •  p  _ 


where  K1  =  150  and  K-,  =  1.75,  and  Um  is  the  superficial  velocity  measured 
at  an  average  pressure  and  based  on  the  mass  flow  entering  the  packed  bed. 
The  coefficients  and  K,  were  determined  by  the  method  of  least  squares, 
from  data  representing  640  experiments  using  spheres,  sand  and  pulverized 
coke  with  air  and  other  gases  including,  C0?,  N2,  CH^  and  gases. 

A  transformation  of  equation  (3)  gives  a  more  conventional  (and  non- 
dimensional)  friction  factor  defined  as. 


F  3  ^ 

v  L 


where  from  equation  (1) 


F  =  150  +  1.75  -7—-, 

v  (.  1  -<t>) 


U  IV 

m  L1 


liquation  (.4)  gives  a  linear  relationship  between  Fy  and  ^  ■ 

In  1951  lirgun  further  refined  this  formula,  equation  (3),  by  con¬ 
ducting  experiments  on  a  single  system  at  various  porosities  [1].  His  re¬ 
sults  included  the  superficial  velocity,  Um  =  4>(Ug  -  U^) ,  which  is  due  to 
the  granular  bed  becoming  mobile  due  to  the  high  pressure,  and  velocity  of 
the  pressure  front.  It  is  based  on  empty  column  cross  section.  His  re¬ 
sulting  formula  was. 

An  1 1  aa-  (U  -  U  )  ,,  ..  (U  -  U  )2d> 

~  =  150  u  — ^ +  1.75  p  -  S  P.  (5) 

<T  Of,  4  b 
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which  can  be  rearranged  as: 

(U  -  U  )  [150(1-$)  *  1.75  -US  P*D.]}  (5a) 


or 


^E  =  P  (U  -  U  )  [150(1-$)  +  1.75  Re]} 


g  P 


(5b) 


Equation  (Sb)  can  be  used  to  define  an  alternate  friction  factor,  i.e.. 


where 


^E  = 


f-  (ug  -  V 

b 


F  =  {  150(1-$)  +  1.75  Re  } 

L  $“ 


(6) 


(6a) 


is  the  non-dimensional  friction  parameter.  The  subscript  'E',  denotes 
the  Ergun  coefficient. 

Equation  (5)  can  be  rearranged  yielding  a  nondimensional  coefficient 
of  drag,  f,  in  which 


_  Ag_  _$ _ 

=  L  g2  (l-$) 


where  G  =  p  UaVg 

Comparing  equation  (7)  with  equations  (6)  and  (6a),  one  sees  that. 


f  = 


Fc 


Re(l-$) 


(8) 


Thus 


b 


(8a) 
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a  correlation  which  is  presented  in  Figure  2,  taken  from  Bird,  Stewart  and 
Lightfoot  [ 4 J .  Here  a  and  b  are  constants  equal  to  150  and  1.75  respecti- 
fully.  This  figure  graphically  shows  Hrgun's  work  blending  the  limits 
of  Blake-Kozeny  equations  and  Burke-Plummer  equations.  The  ordinate  is 
equation  (7)  and  the  abscissa  is  Reynolds  number  defined  in  equation  (4) 
divided  by  the  solids  loading  (l-4>). 

As  stated  above,  Ergun's  correlation  was  tested  to  a  Reynolds  num¬ 
bers  range  of  Re  =  1380.  This  upper  limit  on  the  validity  with  respect 
to  Reynolds  number  limits  the  application  of  this  correlation.  For 

example,  the  flow  of  gases  in  the  ignition  sequence  in  a  ballistic  system 

4  5 

where  the  appropriate  flow  Reynolds  number  may  reach  10  to  10  .  Recently 
Kuo  and  Nydegger  obtained  a  different  pressure  loss  expression  for  gas 
penetrating  the  tightly  packed  beds  of  small  arms  propellant  grains,  at  a 
Reynolds  number,  about  ten  times  as  high  as  those  used  by  Ergun  [5].  They 
attempted  to  simulate  the  Reynolds  number  condition  existing  for  the  igni¬ 
tion  sequence  in  a  gun  cartridge.  Kuo  and  Nydegger  used  a  cold  flow  test 
device  with  a  140  atmosphere  compressor  to  provide  the  high  pressure  gas. 
They  used  actual  propellant  grains  less  than  1mm  in  diameter,  so  that  they 
would  obtain  a  true  packing  configuration.  They  tested  at  Reynolds  numbers 
from  420  to  14,600  and  their  data  provided  the  following  correlation, 

^E=ii41  .  •  [5.05  ♦  276.23]  (9) 

L  0° 


This  expression  can  be  rewritten  as 


L 


(U 

Db  8 


u  ) 

p 


F 


K-N 


(10) 
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where  the  nondimensional  friction  factor,  for  the  Kuo/Nydegger  tests, 
(called  Fj,  )  can  be  expressed  as 

n  2  D  0.87 

fk-n  e  21  [276-23  +  5-05  W  1  (H) 

4> 

Figure  3  compares  F^  and  F^  as  a  function  of  a  Reynolds  number  grouping 
Re 

.  Note  that  Kuo  arrives  at  a  higher  pressure  drop  than  Ergun  at 
low  Reynolds  number,  and  a  lower  drag  than  Ergun  at  high  Reynolds  number. 

In  1978,  Robbins  and  Gough  developed  an  experimental  apparatus  to 
study  the  pressure  drop  through  packed  beds  at  Reynolds  numbers  appro¬ 
priate  to  the  flow  sequence  in  large  caliber  cartridges  [6] .  They  found 

3 

that  at  Re  ^  10  their  measured  friction  factor  for  spheres  agreed  reason¬ 
ably  well  with  Ergun's.  However,  at  Reynolds  number  Re  'v  105,  which  is 
common  to  interior  ballistics,  their  friction  factor  was  57%  of  Ergun's 
extrapolated  results  and  83%  of  Kuo's  extrapolated  values.  They  attempted 
to  provide  agreement  of  their  data  with  the  correlation  of  Kuo  and  Nydegger 
by  utilizing  an  empirical  correction  factor  based  on  the  boundary  layer  de¬ 
veloping  from  the  chamber.  It  is  not  obvious  that  this  logic  is  appro¬ 
priate,  since  in  fact  the  experiments  of  Kuo  and  Nydegger  were  carried  out 
inside  a  very  small  diameter  tube,  while  Robbins  and  Gough  used  a  relatively 
large  diameter  test  section.  Based  on  Benenati  and  Brosilow  [Reference  2] 
and  Figure  lb,  it  seems  logical  that  the  boundary  layer  effect,  as  well 
as  the  best  averaged  uniform  porosity  occurs  when  is  large.  For 

example  Dc/ =  9.4  for  Reference  5,  while  ranged  from  9.6  to  59.9 

in  the  work  reported  in  Reference  6. 

Table  1  summarizes  the  work  done  by  the  above  mentioned  studies. 

Table  2  shows  the  relations  between  the  drag  force  and  the  drag  coeffi¬ 
cient  for  each  study.  Note  that  Robbins  and  Gough  do  not  actually 
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provide  a  correlation,  but  rather  a  method  of  utilizing  measured  pressure 
change  values. 


MOTIVATION 

For  the  past  three  years  work  has  been  underway  at  the  University  of 
Illinois  to  analyze  the  unsteady  reactive  two-phase  flow  associated  with 
DDT  (Def lagration-to-Detonation-Transition)  of  confined  granulated  solid 
propellant  and  explosives.  The  theoretical  models  that  have  been  developed, 
show  there  are  constitutive  relations  which  are  "rate-determining"  func¬ 
tions  [7,8].  In  addition,  these  relations  are  generally  not  well  defined 
for  unsteady  flows  at  high  pressure  and  high  solids  loading,  the  regime 
of  most  interest  in  the  DDT  problem.  One  of  the  most  important  of  these 
constitutive  relations  is  the  gas-particle  drag  interaction.  It  basically 
determines  the  hot-gas  permenability  into  the  unignited  portions  of  the 
granulated  bed.  This  gaseous  flow  resistance  has  a  profound  effect  on 
the  flame  spreading  and  combustion. 

In  order  to  investigate  the  gas-particle  drag  interaction,  a  cold 
flow  experimental  apparatus  was  developed  by  the  author.  In  the  future, 
it  is  planned  to  expand  this  test  facility  into  a  hot  flow,  dynamic  pres¬ 
sure  drop  apparatus. 

The  drag  measurements  are  initially  made  under  steady-state  condi- 

4  5 

tions  at  high  Reynolds  numbers  (10  -  10  ),  and  the  results  will  be  com¬ 
pared  to  work  given  in  References  5  and  6.  The  test  section  was  speci¬ 
fied  to  be  two  inches  in  diameter  compared  to  Robbins  and  Gough's  3  inch 
diameter  section.  The  test  section  size  was  purposely  different  than 
those  utilized  by  Ergun  [1],  Kuo  and  Nydegger  [5],  and  Robbins  and  Gough 
[6],  so  that  one  could  compare  drag  measurements  for  the  same  Reynolds 
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numbers  to  ascertain  the  wall  confinement  effects  which  have  been  briefly 
discussed  above. 

Since  DDT  is  a  transient  phenomenon,  it  will  also  be  desirable  to 
compare  the  steady-state  data  with  transient  pressure  data,  to  determine 
if  the  steady-state  drag  correlations  given  in  Table  2  would  in  fact,  be 
valid  in  a  rapidly  transient  environment.  The  pressure  level  in  the  test 
section  can  be  either  increased  with  time  as  more  gas  is  forced  into  the 
test  section,  or  a  high  pressure  gas  tank  can  be  emptied  through  the  sec¬ 
tion  providing  a  rapid  decay. 

Finally,  if  the  data  at  the  higher  steady  state  Reynolds  numbers  or 
for  the  time-dependent  Reynolds  numbers,  do  not  support  the  current  corre¬ 
lations,  then,  if  possible,  a  new  correlation  will  be  developed. 
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CHAPTER  2 

EXPERIMENTAL  APPARATUS 


Introduction 

This  section  details  the  design  of  the  test  apparatus  and  summa¬ 
rizes  the  test  procedures.  In  designing  the  experimental  apparatus,  there 
were  two  major  limiting  factors,  namely  the  pressure  source  currently 
available,  and  the  funds  available  for  fabrication  and  purchases  of  instru¬ 
mentation. 

The  pressure  source  currently  consists  of  a  bank  of  nine  compressed 
air  tanks  that  supply  a  maximum  working  pressure  of  450  psig,  cold  air 
flow.  The  packed-bed  apparatus  had  to  be  designed  around  this  limitation, 
yet  be  able  to  provide  the  required  mass  flow  and  Reynolds  numbers  of 
interest  to  this  study. 

Figure  4  gives  a  detailed  schematic  of  the  experimental  apparatus 

from  the  plenum  tank  through  the  test  section.  The  eight  inch  long  test 

section,  which  contains  the  packed  bed,  is  two  inches  in  diameter.  As 

stated  previously,  Kuo's  and  Nydegger's  results  were  obtained  from  a  0.30 

inch  diameter  test  section,  D  ,  allowing  the  ratio  D  /D,  =  9.4.  Accord- 

c  CD 

ing  to  Benenatic  and  Brosilow  [2],  it  seems,  that  appropriate  average 
porosity  and  minimum  boundary  layer  effect  occur  with  large  Dc/DjJ  values, 
used  for  example  by  Robbins  and  Gough  [6].  Hence,  a  test  section  diameter 
larger  than  Kuo's  and  Nydegger's  [5],  but  smaller  than  Robbin's  and  Gough's 
was  designed,  since  this  would  enable  us  to  compare  the  drag  measurements 
for  identical  conditions  of  particle  size,  mass  flow  and  hence  Reynolds 
number,  but  for  different  ratios  of  D^/D^ . 


The  eight  inch  length  of  the  test  section  was  similar  to  other 
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researchers,  and  allows  enough  distance  from  the  ends  of  the  test  section 
to  the  first  and  last  pressure  transducer,  so  that  the  entrance  and  exit 
flow  patterns  do  not  bias  the  pressure  data.  Also  the  length  chosen 
would  diminish  any  errors  that  may  occur  due  to  nonuniformity  in  packing. 

The  two  inch  distance  between  pressure  transducers  was  arbitrarily  chosen, 
but  did  allow  for  appropriate  mounting  of  the  gauges  on  the  outside  of  the 
pipe  test  section.  Three  pressure  transducers  were  chosen  so  that  one  could 
record  at  least  two  to  three  pressure  gradients  in  the  flow. 

Details  of  Specific  Hardware  Utilized 

Figure  5  is  a  schematic  of  the  entire  experimental  apparatus.  Bour¬ 
don  gauge  #1  is  a  Solfrunt  4.5  inch  test  gauge  with  a  range  of  0-1000  psig 
and  a  +_  0.25%  of  span  accuracy,  its  purpose  is  to  continuously  monitor 
the  supply  pressure.  This  is  a  moderately  priced  gauge,  yet  it  has  the 
range  to  be  used  with  a  higher  pressure  source.  A  Grove  Pressure  Reduc¬ 
ing  Regulator,  Model  82-829  (720  psig  max.  inlet),  allows  us  to  regulate 
the  compressed  air  source  to  the  desired  test  pressure.  This  regulator 
could  handle  the  source  pressure,  required  mass  flow,  and  in  addition, 
regulate  itself  so  that  the  desired  test  pressure  in  the  plenum  tank 
would  not  fluctuate.  The  1-1/4  NPT  pipe  was  compatable  with  the  Regu¬ 
lator  and  mass  flow  requirement. 

The  regulated  compressed  air  enters  a  plenum  tank  with  an  inner 

diameter  of  2.99  inches.  This  enlarged  diameter,  from  1.27  inch  in  the 

pipe,  to  2.99  inch  in  the  plenum,  was  designed  to  bring  the  gas  to  a 

very  low  velocity.  The  three  straightening  vanes  are  to  smooth  the  flow 

of  air  in  the  plenum.  Bourdon  gauge  #2  is  the  same  type  as  Bourdon  gauge 

#1.  Gauge  #2  indicates  the  test  pressure,  equal  to  the  stagnation  pressure 

P  ,  which  is  used  in  the  calculation  of  mass  flow.  At  the  end  of  the 
o 
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plenum  tank  is  a  Sonic  Flow  Nozzle  from  the  American  Meter  Division,  Singer 
Corporation.  This  nozzle  has  a  throat  diameter  of  0.375  inches,  and  an 
accuracy  of  +_  0.15%  at  10  psig.  This  nozzle  was  selected  because  of  its 
successful  use  by  Robbins  and  Gough  [6],  its  moderate  cost,  and  its  opera¬ 
tion  to  provide  the  required  mass  flow  rates.  This  nozzle,  a  converging- 
diverging  nozzle,  provided  supersonic  flow  at  its  exit  as  explained  in 
Appendix  B.  The  high  speed  flow  is  then  gently  decelerated  so  that  lami¬ 
nar  flow  enters  the  test  section. 

Opposite  to  the  Bourdon  gauge  #2,  is  an  Omega  Thermocouple.  It  is  a 
chromel- alumel ,  CASS-18G-12,  grounded  junction  probe,  that  is  designed  to 
provide  fast  response  under  high  pressures.  This  thermocouple  measures  the 
stagnation  temperature,  Tq  (in  the  plenum  tank)  required  for  use  in  the 
mass  flow  calculation. 

The  sonic  nozzle  opens  into  an  expansion  section  which  connects  into 
a  straightening  vane  section.  Based  on  the  recommended  minimum  lengths  of 
pipe  preceding  and  following  orifices,  flow  nozzles  and  venturies  from  the 
ASME  "Fluid  Meters,"  Chart  [9],  a  nine  inch  long  expansion-straightening 
section  would  be  required.  After  running  several  tests  in  this  configura¬ 
tion,  it  was  decided  to  lengthen  the  expansion  section,  so  that  the  flow 
would  definitely  be  fully  developed  when  it  entered  the  test  section. 

All  straightening  vanes  are,  in  fact,  grids,  constructed  of  1/4  inch 
ground  stock,  so  that  the  straightening  vane  is  a  mesh  grid  only  1/4  inch 
wide.  Figure  6  is  a  photograph  of  one  of  the  mesh  grids  as  inserted  into 
the  straightening  section.  Note  that  several  sets  are  positioned  in  the 
respective  sections.  The  straightening  vanes  were  constructed  in  this  man¬ 


ner  so  that  boundary  layer  interaction  would  be  minimized. 

The  eight  inch  long  and  two  inch  diameter  test  section  is  directly 
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mounted  to  the  straightening  section.  This  section  has  a  screen  mesh  at 
the  entrance  and  a  metal  stainless  steel  grid  at  the  exit,  to  secure  the 
packed  bed  in  the  test  section  as  shown  in  Figure  7.  There  are  three  rapid 
response  pressure  transducers  mounted  on  the  test  section  to  measure  the 
pressure  drop  of  the  gas  as  it  flows  through  the  packed  bed.  The  trans¬ 
ducers  are  Setra  Systems,  Model  205,  with  a  pressure  range  0-500  psig,  a 
one  millisecond  response  time  and  an  accuracy  of  0.11%  at  full  scale.  This 
type  transducer  provides  output  signals  for  detection  of  any  fluctuation 
in  the  steady-state  air  flow,  and  for  continuous  monitoring  of  the  pres¬ 
sure  in  the  transient  tests. 

The  Omega  thermocouples,  as  described  previously,  monitor  the  tempera¬ 
ture  in  the  test  section  and  also  after  the  gas  exits  the  test  section. 

All  three  thermocouples  are  connected  to  an  Omega  Digital  Thermometer, 

Model  2160A,  with  a  response  time  of  less  than  2  seconds,  and  a  maximum 
error  including  NBS  conformity  of  +_  °1F.  The  Omega  Thermocouples  and  digi¬ 
tal  thermometer  were  on  hand  and  were  used  because  they  could  accurately 
and  rapidly  record  the  temperatures. 

For  the  decay  transient  experiments  the  regulator  is  not  used.  In¬ 
stead  a  high  pressure  flex  line  (7000  psig  limit)  is  connected  one  inch 
downstream  from  the  regulator  onto  the  pipe,  as  shown  in  Figure  8.  This 
line  is  connected  to  a  high  pressure  two  way  ball  valve  (3000  psig  limit), 
which  is  connected  to  a  compressed  air  bottle  (2400  psig).  Bourdon  gauge 
K2  is  replaced  by  a  Setra  System,  Model  204E  pressure  transducer,  (0-1000 
psig) ,  to  continuously  monitor  the  pressure  in  the  plenum  tank.  The  thermo¬ 
couple  in  the  plenum  tank  is  connected  to  a  storage  oscilloscope,  as  de¬ 
scribed  below.  This  allows  the  continuous  monitoring  of  the  temperature 
in  the  plenum  tank,  as  the  pressure  decays  or  increases,  whichever  the 
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case  may  be.  Figure  9  and  10  show  the  exterior  of  the  apparatus  ready  for 
use  in  the  steady  state  tests.  Figure  11  shows  the  arrangement  for  the 
transient  experiments. 

The  data  from  the  Setra  Transducers  is  recorded  on  a  Tektronic  5100 
series  storage  oscilloscope,  and  the  traces  are  permanently  recorded  with 
a  Tektronix  C-S  oscilloscope  camera.  Thermo  Systems,  Model  1076,  digital 
RMS  voltmeters,  were  also  used  as  a  back  up  and  comparison  check  when  re¬ 
cording  the  steady  state  output  from  the  transducers. 

The  wall  thickness  of  the  various  pipes  and  tubing  were  determined  by 
performing  a  hoop  stress  analysis  of  each  individual  section,  so  that  each 
section  could  withstand  a  maximum  pressure  of  1000  psi  with  a  safety  fac¬ 
tor  of  five. 

Reduced  Chamber  Diameter 

In  order  to  determine  the  effect  on  the  measured  pressure  gradient  of 
the  wall  boundary  layer  build  up  and  possible  channeling  of  flow  through 
the  "annular"  region  at  the  wall,  a  sleeve  was  designed  to  be  inserted  into 
the  test  section.  The  metal  sleeve  was  press  fit  into  the  section,  such 
that  the  pressure  port  holes  were  aligned  with  the  original  pipe  section. 

The  insert  was  one  inch  in  diameter,  reducing  the  flow  area  of  the  standard 
two  inch  diameter  test  section  by  one  fourth.  This  was  done  so  that  measure¬ 
ments  could  be  carried  out  at  two  different  ratios  of  D  /D,  ,  for  the  same 

c  b 

size  beads.  Discussions  were  given  in  the  previous  chapter  which  indicate 
why  this  type  of  data  may  be  required.  Table  3  summarizes  the  components 
for  the  packed  bed  test  facilities  used  by  Kuo  and  Nydegger  [5],  Robbins 
and  Gough  [6],  and  this  study. 
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Calibration  and  Experimental  Procedure 


The  Bourdon  gauges  and  pressure  transducers  were  all  calibrated  on 
a  dead  weight  tester  prior  to  being  installed  on  the  experimental  appara¬ 
tus.  These  calibrations  compared  favorably  with  those  supplied  by  the 
gauge  manufacturer.  See  Appendix  C  for  calibration  details. 

The  experimental  procedures  for  conducting  steady  state  and  tran¬ 
sient  experiments  can  be  found  in  Appendix  D.  The  most  important  factor 
when  conducting  the  experiments  was  to  compact  the  test  bed  at  a  higher 
pressure  than  would  be  used  during  the  tests.  This  ensures  that  the  packed 
bed  will  not  further  compact  during  the  tests. 

The  experimental  data  obtainable  from  this  apparatus  is  limited  by 
the  compressed  air  source  of  450  psig,  unless  high  pressure  air  bottles 
are  manifolded  into  the  plenum.  In  the  latter  case  the  time  for  carrying 
out  steady-state  (high-pressure)  experiments  is  limited  as  these  tanks  empty. 

The  input  pressures  available  limit  the  flow  rate  and  consequently  the 
Reynolds  number  range  that  presently  can  be  studied.  It  is  not  now  possi¬ 
ble  to  obtain  heated,  high  pressure  air  with  this  apparatus. 

Specific  Experiments  Performed 

Flow  resistance  experiments  were  performed  mainly  using  the  two  inch 
diameter  test  section  in  which  spherical  beads  of  0.96,  3.0,  and  6.0  mm 
diameter  were  packed.  The  pressure  gradient  data  was  reduced,  analyzed, 
and  compared  to  the  theoretical  correlation  of  Ergun,  Kuo  and  Nydegger,  as 
well  as  to  the  experimental  data  of  Robbins  and  Gough. 

Table  4  summarizes  almost  all  of  the  steady-state  experiments  per¬ 
formed.  The  combination  mixture  of  6  mm  and  0.96  mm  beads  was  made  to  de¬ 
termine  any  correlation  between  pressure  drop,  porosity  and  size  mixing. 


State  Experiments 


Cylinders  100,  150,  200, 

L=0.403  in.  47g  0>443  250,300,350,  16595  <  Re  <  59251 

D  =.183.  400 
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(.Obviously,  many  more  mixture  ratios,  small/large,  could  be  tested.) 

The  packed  bed  of  inert  cylindrical  propellant  was  tested  at  the 
same  pressures  used  for  the  spherical  beads.  The  pressure  drop  data  was 
used  to  verify  the  observations  made  by  Robbins  and  Gough  [6]  about  the 
geometrical  scaling  of  particles  on  the  bed  resistance. 

The  following  chapter  summarizes  the  data  and  the  conclusions  that 
can  be  drawn  from  these  steady -state  experiments. 

Transient  Pressure  Input 

There  are  two  type  transients,  progressive  and  decay.  The  progres¬ 
sive  transient  occurs  when  the  pressure  in  the  test  section  is  increased 
from  0  psig  to  400  psig  very  rapidly,  by  opening  the  regulator.  The  de¬ 
cay  transient  occurs  when  a  high  pressure  air  bottle,  which  is  attached 
to  the  apparatus  as  previously  explained,  is  allowed  to  empty  from  2400 
psig  to  the  room  atmospheric  pressure. 

The  time  for  the  progressive  transient  depends  only  on  how  rapidly 
one  is  able  to  completely  open  the  regulator  by  hand.  Presently,  it  takes 
from  one  to  two  seconds. 

The  time  for  the  decay  transient  depends  on  two  factors.  First,  the 

total  volume  of  the  air  in  the  bottle  will  dictate  how  long  it  takes  to 

flow  through  an  opening.  Typical  commercially  available  bottles  of  high 

3  3 

pressure  air,  are  tanks  with  volumes  of  309  ft  or  296  ft  .  Secondly, 
the  size  of  the  opening  in  the  air  bottle  will  determine  how  much  air  can 
pass  out  of  the  bottle  per  second.  In  our  specific  case  the  manufacturer 
of  the  air  bottle  has  specified  the  size  of  the  opening  from  the  bottle, 
and  this  can  not  be  changed. 


In  Appendix  E,  the  method  for  calculating  the  time  to  empty  a  pres¬ 
sure  chamber  is  explained.  It  takes  2.73  minutes  (164  seconds)  for  our 
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decay  transient,  using  an  air  bottle  with  309  ft3  of  air  at  an  initial 
pressure  of  2400  psig  and  at  room  temperature. 

The  transient  pressure  output  at  each  of  the  three  pressure  trans¬ 
ducers  in  the  test  section  were  recorded  on  an  oscilloscope.  This  con¬ 
tinuous  pressure  gradient  data  for  the  0.96  mm,  3  mm  and  6  mm  diameter 
beads  was  compared  to  the  discrete  steady-state  data,  to  determine  whether 
the  drag  in  the  steady -state  case  was  the  same  as  that  in  the  transient  case. 
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CHAPTER  3 


ANALYSIS  OF  TEST  RESULTS 


For  the- test  conditions  as  outlined  in  Table  4  (Chapter  2)  this  coeffi¬ 
cient  was  calculated  utilizing  the  above  relation,  where  U  is  deter- 

avg 

mined  by  the  mass  flow.  (Refer  to  Appendix  B.)  Table  2  lists  two  specific 
formulas  for  the  drag  coefficient  as  reported  in  the  literature  by  Ergun 
[3]  and  Kuo  and  Nydegger  [5]. 

Figure  12  plots  these  correlations  as  a  function  of  Reynolds  number 
as  lines  of  constant  porosity,  4>.  It  should  be  noted  that  the  correlation 
by  Ergun  has  been  extrapolated  well  beyond  the  range  of  Reynolds  number 
where  this  correlation  was  developed.  For  the  two  inch  interior  diameter 
test  section,  the  calculated  F  (from  equation  (12)]  is  plotted  for  three 
initial  size  spherical  beads,  The  porosities  for  each  test  grouping  is 
noted  on  the  graph. 

It  would  appear  that  for  the  larger  particles  (3mm  and  6mm),  the 
correlation  by  Kuo  and  Nydegger  is  adequate.  Unfortunately,  the  sub-milli¬ 
meter  beads,  provided  data  that  gives  Fv  significantly  less  than  that  corre¬ 
lation,  at  least  in  the  range  of  Reynolds  number  tested  here.  This  is 
surprising  since  the  work  reported  in  Reference  5  utilizes  sub-millimeter 


diameter  beads,  (0.826mm).  As  was  shown  in  Table  3  of  the  previous  chapter. 
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the  ratio  of  test  chamber  diameter  to  bead  diameter,  D^D^,  for  the  Kuo 
and  Nydegger  experiments  was  a  value  less  than  ten.  For  the  0.96mm  spheri¬ 
cal  beads  used  here,  this  ratio  was  greater  than  fifty.  It  should  be 
noted  that  when  using  the  6mm  diameter  beads,  the  ratio  l^/D^  is  eight 
and  a  half,  a  number  approaching  that  used  in  Ref.  5.  Note  that  this  con¬ 
dition  correlates  well  with  the  published  formula  for  F  . 

In  order  to  more  clearly  determine  the  effect  of  the  containers,  i.e., 
the  test  section  diameter,  a  smaller  test  section,  exactly  one  inch  in¬ 
terior  diameter  was  utilized.  Table  4b  has  summarized  the  tests  performed 
with  that  modified  apparatus.  Since  the  area  of  this  section  is  now  one 
fourth  that  of  the  two  inch  test  section,  one  achieves  much  higher  average 
velocities  for  the  same  mass  flows.  Therefore,  the  range  of  the  test 
Reynolds  number  is  greatly  extended.  The  Reynolds  number  range  now  avail¬ 
able  was  well  above  those  in  the  tests  in  Reference  5.  Figure  13  presents 
the  coefficient  of  drag,  F  versus  Reynolds  number  for  the  smaller  diam- 
eter  test  section.  The  data  clearly  indicates  that  the  correlation  pro¬ 
posed  by  Kuo  and  Nydegger  is  no  longer  valid  for  this  flow  regime.  It  now 
appears  that  there  are  no  particles  within  the  range  tested,  that  would 
fit  the  correlation  of  Reference  5,  as  there  were,  when  the  larger  diam¬ 
eter  test  section  was  utilized. 

Robbins  and  Gough  [oj  had  also  reported  that  their  pressure  drop  data 
could  only  be  marginally  described  by  the  correlation  of  Kuo  and  Nydegger. 
Their  data  also  indicates,  as  was  done  here,  that  the  correlation  of  Ergun 
should  not  be  used  at  these  large  values  of  Reynolds  numbers.  Some  of 
the  data  reported  by  Robbins  and  Gough  will  be  compared  with  our  findings 
in  a  following  section. 
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Search  for  High  Reynolds  Number  Drag  Coefficient 

When  reviewing  the  type  of  data  shown  in  the  previous  two  figures, 
it  becomes  evident  that  one  does  not  arrive  at  one  unique  coefficient  of 
drag  at  a  fixed  Reynolds  number  when  the  test  sections  provide  for  a  dif¬ 
ferent  flow  area.  That  is,  for  any  one  fixed  size  spherical  bead,  a  plot 
of  Fv  versus  Reynolds  number  gives  a  relation  that  has  a  significantly 
different  slope.  Data  shown  on  Figure  14  for  the  0.96mm  bead  [Part  (a)] 
and  the  3mm  bead  [Part  (b)],  clearly  bear  this  out. 

The  variation  shown  in  Figure  14  is  due  to  the  fact,  that  at  the  same 
Reynolds  number,  the  average  velocity  is  significantly  different.  Equal 
Reynolds  number  at  different  average  velocities  is  achieved  by  adjusting 
the  average  gas  density  (or  pressure).  Since  Re  =  m  D^/A  0*  this  is 
achieved  by  altering  the  mass  flow  to  match  the  change  in  pipe  area.  One 
must  conclude  that  Fy  is  not  solely  a  function  of  the  flow  Reynolds  num¬ 
ber  as  has  been  summarized  by  previous  investigators  and  as  reviewed  in 
Chapter  1. 

After  inspection  of  the  pressure  drop  data,  it  was  found  that  for 

any  given  size  bead,  irrespective  of  the  test  section  flow  area,  Fy  corre- 

2 

lates  as  a  function  of  Re/Uavg  .  This  is  shown  in  Figure  15,  indicating 

that  the  data  given  in  Figure  14  collapses  to  one  line  when  Fy  is  plotted 

2 

as  a  function  of  Re/Uavg  .  This  correlation  is  only  approximately  true 
for  the  tests  with  the  6mm  spheres,  but  this  is  probably  due  to  the  fact 
that  for  the  one  inch  interior  diameter  test  section,  the  6mm  beads  are 
too  large  to  represent  a  uniformly  packed  bed. 

Further  replotting  of  the  data  showed  that  the  effect  of  the  bead 

2 

size  can  also  be  correlated  if  one  plots  F  as  a  function  of  bD,  Re/U  , 

r  v  b  avg  ’ 

where  b  =  a  constant.  That  correlation  is  shown  in  Figure  16.  Here  the 
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Figure  16.  Coefficient  of  Drag  as  a  Function  of  Bead  Size 
Reynolds  Number  and  Velocity. 


6mm  bead  data  (with  the  2  in.  I.D.  test  section)  also  correlates  with  this 
functional  relation. 

This  surprising  correlation  of  coefficient  of  drag  with  the  flow  para¬ 
meter  indicates  that  the  flow  velocity,  or  more  correctly,  the  kinetic 

2 

energy  of  the  flow  (n,  )  is  a  rate  determining  factor  for  Fy.  Figure 

17  plots  Fy  versus  Reynolds  number  only.  The  data  shown  here  is  identical 
to  that  given  in  Figure  14b,  except  that  the  data  of  Robbins  and  Gough 
for  a  sphere  nearest  the  3mm  bead  tested  here,  is  added.  As  has  been  in¬ 
dicated  in  Table  3  (Chapter  2),  Robbins  and  Gough  carried  out  their  experi¬ 
ments  in  a  three  inch  interior  diameter  pipe.  Table  5  below  summarizes 
the  range  of  average  velocities  tested  with  the  sections  used  here  and  with 
the  apparatus  used  by  Robbins  and  Gough,  and  Kuo  and  Nydegger. 

Table  5.  Average  Flow  Velocities 


Investigator 

Test  Section 
Diameter  (inch) 

Range  of 
U,vg  <FI*> 

Wilcox/ Krier 

1.00 

305  -  340 

Wilcox/ Krier 

2.00 

83  -  132 

Robbins/Gough 

3.00 

54  -  98^ 

Kuo/Nydegger 

0.3059 

23  -  121^ 

Estimated 


It  appears  that  for  moderately  low  velocities  (below  100  fps) ,  a  unique 
correlation  between  and  Reynolds  number  is  possible.  But  for  the 
higher  average  gas  velocity,  a  correlation  of  the  nature  indicated  in 
Figure  16  is  more  likely. 

An  alternate  way  to  show  that  Reynolds  number  alone  does  not  corre¬ 


late  the  resistance  to  flow  in  a  packed  bed,  is  to  plot  an  alternate 
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Reynolds  Number  (xIO'3) 


Figure  17.  Comparison  of  Drag  Coefficient  Data,  Indicating  Flow 
Velocity  Effect. 
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coefficient  of  drag,  f,  as  defined  in  Hqs.  (7)  and  (8).  The  data  pre¬ 
sented  by  Robbins  and  Gough  is  represented  in  this  form.  Clearly  a  ve¬ 
locity  effect  other  than  that  given  in  the  Reynolds  number,  dominates 
this  altered  form  of  the  coefficient  of  drag,  as  shown  in  Figure  18,  where 
£  is  plotted  versus  Reynolds  number. 

Formulation  of  High  Reynolds  Number  Drag  Coefficient 

The  classical  equation  for  the  coefficient  of  drag  is  =  a  +  b  Re, 
based  on  Ergun's  work  [1],  or  more  specifically  Eq.  (4),  Fy  =  150  +  1.75* 
Re/(1-<J>).  As  discussed,  the  data  at  the  higher  Reynolds  numbers  indi¬ 
cate  that  this  type  of  linear  relation  with  Reynolds  number  is  not  cor¬ 
rect.  Instead,  one  notes  that  for  the  current  study, 

r  .Re 


U 


(13) 


avg 


Upon  comparing  Eq.  (13)  with  the  data,  it  was  found  that  b^  was  a  function 
of  bead  diameter.  A  plot  of  b  versus  on  log- log  paper  resulted  in  a 
fairly  accurate  condition  that 


b  =  a  D, 


where  a  =  3.33  x  106  ft/sec2  •  Equation  (13)  then  becomes 


(14) 


(15) 


where  Fy  is  nondimensional ,  as  required.  Recall  that  Eq.  (4a)  defines 
the  Reynolds  number,  and  Eq.  (B-3)  the  average  gas  velocity/  uaVg»  as 


Re  =  p 


U  0  D. 
avg  b 


avg 


P  <(>  A  . 

pipe 


Substituting  Eq.  (15)  into  Eq.  (12),  the  pressure  gradient  becomes. 


^  This  correlation  was  specifically  developed  at  moderate  to  high  gas  ve¬ 
locities,  and  for  Reynolds  number  greater  than  10^.  It  should  not  be 
extrapolated  to  very  low  velocities,  since  Fy-*  00  as  Ua  ■+•  0. 


obis—  ro _ ot  »  »  -i  oo  isO _ ro _ o'  *■  ai  o>  >igioO 
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Wilcox  /  Krier  -3  mm  Beod 
A  =  7 7  in2  Pipe  Data 
A  =  ^4  in2  Pipe  Data 


Robbins /  Gough  -  4.75  mm  Bead 
•  =  %  7 T  in2  Pipe  Data 


Kuo  /  Nydegger 


.01 


Wilcox /Krier 
High  Velocity  {  *- 
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Figure  18.  Alternate  Form  of  Coefficient  of  Drag,  Indicating  Velocity  F.ffect. 


45 


Eg.  f 163  gives  AP/L  in  units  of  lbf/ft' 


6  2 

Recall  the  constant  a_  equals  3.35  x  10  ft/sec  . 


(16) 


In  SI  units  Eq.  (16)  is  simply  expressed  as 


—  =  a'  p  <J>  (16a) 

3  3 

where  AP/L  is  in  units  of  nt/m  ,  p  in  units  of  kg/m  ,  and  a'  equals 

6  2  3 

1.015  x  10  m/sec  .  Again  Eq.  (16)  should  only  be  utilized  for  Re  >  10  . 

Equation  (16)  was  derived  from  data  utilizing  almost  exclusively 
uni-sized  spherical  particles  packed  randomly  to  provide  the  highest  solid 
loadings  (or  lowest  porosity,  <j>) .  The  porosity  range  was  therefore  limited 
to  0.38  <  $  <0.43.  As  such,  the  relation  [Eq.  (16)]  should  not  be  used  to 
indicated  the  sensitivity  of  the  flow  resistance  as  a  function  of  solids 
loading. 

As  a  matter  of  fact,  when  a  bi-modal  mixture  of  6mm  and  0.96mm  beads 
was  mixed,  about  50/50  by  weight,  the  resulting  porosity  was  <p  *  0.29. 

For  the  same  mass  flow  through  the  bed,  (as  for  either  the  all  6mm  or  all 
0.96mm  beads),  the  data  clearly  indicated  that  the  pressure  gradient,  AP/AL 
was  larger  for  the  bi-modal  mix,  since  the  density  of  the  gas  was  addition¬ 
ally  larger. 

That  the  correlation  given  by  Eq.  (1C)  is  fairly  accurate  can  be  seen 
by  reviewing  the  data  for  AP/AL  given  in  Figure  19,  where  the  equation  is 
compared  to  the  measure  pressure  drop. 

It  should  be  noted  that  there  is  some  similarity  between  our  Eq.  (16) 
and  the  standard  fluidized  bed  equation.  That  is,  a  vertical  bed  becomes 
fluidized  when  a  gas  of  lower  density,  than  that  of  the  bed  particles. 


J  0  10  20  30  40  50  60  70 

^  calculated  (xIO'3)  ,  Eqn.[|(>] 


Ficure  19.  Comnarison  of  Measured  Pressure  Drot>  to  Predicted  Pressure  Drcn. 
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flows  upward  through  the  bed  causing  the  particles  to  no  longer  rest  on 
each  other,  but  be  suspended  by  the  upward  force  of  the  gas.  As  the  gas 
velocity  increases,  the  particles  will  expand.  This  fluidized  state  is  ex¬ 
pressed  by  the  following  equation  [11], 

=  g  [4>po  +  U-4>)Pb]  (17) 

The  horizontal  packed  bed  in  this  work  did  not  become  fluidized  be¬ 
cause  it  was  tightly  packed  and  then  the  bed  was  retained  by  metal  disks 
at  each  end  of  the  test  section.  Much  more  analysis  would  be  necessary 
in  order  to  explain  the  correlation  obtained  at  high  velocities  [Eq.(16)] 
to  a  similar  correlation  for  fluidized  beds  [Eq.  (17)]. 

As  one  can  see  from  Eq.  (16),  the  pressure  gradient  is  proportional 
to  the  average  gas  density  in  the  test  section.  Rearranging  Eq.  (16), 
yields 

AP/Ax  =  =  cp  (18) 

8c 

where  Cj  =  a<t>/gc  RT,  which  is  approximately  constant  since  the  gas  tem¬ 
perature,  Tj,  varies  only  by  one  or  two  degrees  in  the  region  where  the 
pressure  is  measured.  Integrating  Eq.  (18)  gives 

£n  P  =  C^x  (19) 

Figure  20  is  a  plot  of  Eq.  (19),  and  it  can  be  seen  that  the  slope,  C^, 
is  constant  independent  of  initial  gas  pressure.  In  fact,  the  slope  will 
be  constant  irregardless  of  the  bead  diameter  and  test  section  used.  Figure 
21  graphically  illustrates  this  point.  A  plot  of  Robbin's  and  Gough's  data 
in  the  form  of  Eq.  (19)  yields  the  same  results,  as  shown  in  Figure  22 
[12]. 

Further  Analysis  of  Developed  Correlations 


From  Bird,  Stewart  and  L’ghtfoot  [4]  for  steady  flow  where  Dary's 


Apipe  =  77/^  in2,  A  =  3  mm 


Law  for  gas  permeability  is  assumed  valid,  (low  Reynolds  number  range), 
one  can  show  that  for  steady  isothermal  gas  flow 


V2o2  =  0  or  V2P2  =  0  (20) 

which  for  one  dimensional  flow  implies 

P  ~  /“x  (21) 

To  determine  how  P  should  vary  with  distance,  for  the  moderate  Reynolds 
number  flows,  where  the  coefficient  of  friction 


F  =  a  +  b  Re  or  F  ~  C_  Re 
v  v  3 


(22) 


one  can  substitute  Eq.  (22)  into  (12)  to  get 


f  =  c  Re  ^8 


(23) 


and  since. 


m 


avg  p  A  . 
s  pipe 


and  Re  = 


A  . 
pipe 


then  simplifying  Eq.  (23)  yields 


dP.  r  1 

dx  L4  P 


(23a) 


2  2 

where  =  C^m  RT/Ap..  <t>  D^,  and  integrating  Eq.  (23a)  gives 


P  =  C5X  or 


^7 


(23b) 


as  indicated  above  in  Eq.  (21).  The  expression,  Eq.  (23a),  assumes  mass 
flow  (m)  is  an  independent  parameter.  If  instead  one  assumes  that  m  ~  P, 
which  was  the  case  in  this  work,  then  Eq.  (23)  becomes,  after  integration, 


in  P  =  C,X 
o 


(24) 


which  is  the  same  form  obtained  in  Eq.  (19).  In  conclusion,  the  pressure 
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gradient,  as  expressed  in  Eq.  (16)  and  the  coefficient  of  drag  as  expressed 
in  Eq.  (15)  are  logical  formulations  for  these  expressions  at  high  Reynolds 
numbers . 


Transient  Experiments 

As  discussed  in  Chapter  2  the  transient  tests  were  conducted  by  empty¬ 
ing  high  pressure  air  into  the  test  section  from  a  standard  commercial  tank. 
It  was  found  that  the  time  for  this  decay  transient  was  2.70  to  3.0  minutes. 
Appendix  E  shows  why  this  was  so.  The  time  period  for  a  progressive  tran¬ 
sient  experiment  ranged  from  1.5  to  5.3  seconds.  Both  the  decay  and  pro¬ 
gressive  transients  were  much  longer  than  the  desired  millisecond  time 


transients. 


The  transient  tests  were  conducted  using  only  the  two  inch  test  sec¬ 
tion.  The  resulting  flow  resistance  for  the  progressive  transient  tests 
varied  only  slightly  from  the  steady  state  tests,  as  shown  in  Figure  23. 
With  the  available  data  one  can  conclude  that  the  derived  formula  for  the 


coefficient  of  drag,  Eq.  (15),  and  pressure  drop,  Eq.  (16),  will  accurately 
predict  the  transient  results.  This  can  be  explained  by  comparing  the  total 
flow  transient  period,  which  ranged  from  1.5  to  5.3  seconds,  with  a  charac¬ 
teristic  "time",  Tl,  defined  as,  T^  =  Ax/Uavg.  The  distance  Ax  is  1/6  foot 
(pressure  gauge  separation),  and  with  ranging  from  80  fps  to  122  fps, 

r  ranges  from  2  milliseconds  to  1.36  milliseconds.  Unless  the  flow  experi¬ 
ment  could  be  carried  out  to  change  the  pressures  significantly  in  time 
periods  of  the  order  of  T^,  one  should  expect  that  the  correlations  de¬ 
veloped  in  these  (steady-state)  experiments  would  be  valid  for  moderately 
transient  conditions. 


Additional  Tests  with  Various  Size  Particles 

A  mixture  of  6mm  beads  and  0.96mm  beads  was  tested.  The  mixture 
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contained  50.6%,  6mm  beads,  and  49.4%,  0.96mm  beads,  by  weight.  The  effec¬ 
tive  diameter^  of  the  mixture  particles  was  1.249mm  (0.049  in.).  The 
porosity  of  this  mixture  was  0.2962  which  indicates  that  the  0.96mm  beads 
filled  the  voids  between  the  6mm  beads  causing  a  very  tight  packing  of  the 
particles.  A  plot  of  versus  Reynolds  number  for  this  mixture  gave  a 
line  which  was  slightly  above  that  for  the  0.96mm  bead,  which  one  would  ex¬ 
pect  since  the  effective  diameter  of  the  particles  was  larger. 

Ottawa  Sand  was  also  packed  into  the  two  inch  test  section.  The  effec¬ 
tive  diameter  of  the  sand  was  0.7315mm  (0.0287  in.).  A  plot  of  the  result¬ 
ing  Fy  versus  Reynolds  number,  gave  a  line  slightly  below  that  for  the 
0.96mm  bead,  which  was  shown  in  Figure  14a. 

The  results  of  the  mixture  and  sand  tests,  indicate  that  the  0.96mm 
bead  data  is  accurate.  Also,  the  equations  for  F^  and  pressure  gradient, 
(15)  and  (16)  respectively,  would  correctly  predict  the  results  of  the 
above  tests. 

Finally,  a  test  was  carried  out  in  the  two  inch  test  section  using  a 
packed  bed  of  cylindrical  shaped  solids,  10.23mm  (0.4028  in.)  long  and 
4.65mm  (0.183  in.)  in  diameter.  The  resulting  values  of  the  coefficient 
of  drag,  F^  when  plotted  versus  Reynolds  number,  are  in  agreement  with 
the  results  Robbins  and  Gough  [6]  obtained,  if  an  effective  diameter  of 
=  0.183  inch  is  used. 


4  (effective) 


=  2 


r(1-*bed) 

Lo-*,) 
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N. 

l 


N, 


Total 


1/3 


where  N  =  number  of  particles. 


Future  Work  Required 

The  formulated  correlation  to  predict  pressure  drop,  Eq.  (16) 

[AP/L  =  a  p  <}>/ gcJ,  was  compared  to  some  actual  pressure  gradient  data 
from  Robbin's  and  Gough's  work  [12].  It  was  found  that  the  value  of  a 
for  their  data  was  2.56  x  10  ft/sec  instead  of  3.33  x  10°  ft/ sec  as 
developed  in  this  work.  Hence,  the  predicted  pressure  gradients  were 
23%  higher  than  the  measured  pressure  gradients  at  their  high  Reynolds 

3 

number  range  (Re  '  80  x  10  ). 

In  light  of  these  results,  further  work  should  be  conducted  in  order 
to  test  the  validity  of  the  correlations  for  the  coefficient  of  drag,  F.q. 
(IS),  and  the  pressure  gradient,  Eq.  (16).  This  would  entail  construct¬ 
ing  an  additional  three  inch  test  section  and  obtaining  a  pressure  reser¬ 
voir  which  would  increase  the  inlet  pressure  to  the  packed  bed  by  at 
least  an  order  of  magnitude.  To  fabricate  the  high  pressure  vessel  of  a 
significant  capacity,  here  at  the  University  of  Illinois,  as  well  as  the 
purchase  of  a  compressor  that  could  replenish  the  reservoir  in  the  order 
of  several  thousand  pounds  per  square  inch,  would  require  a  substantial 
investment. 

An  alternate  method  of  increasing  the  pressure  in  the  test  section, 
is  to  manifold  several  commercial  gas  bottles  and  reduce  the  nozzle’s 
throat  area,  resulting  in  a  high  stagnation  pressure  in  the  plenum  tank 
and  high  pressures  in  the  test  section.  Tests  would  be  conducted  over 
the  same  Reynolds  number  range,  but  with  higher  pressures.  For  safety 
reasons  this  was  not  attempted. 

Conclusions 


The  work  reported  here  indicates  that  in  gaseous  flow  through  packed 
beds,  Ergun's  [1]  correlation  is  valid  for  Reynolds  number  less  than  1400. 
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As  the  Reynolds  number  increases,  there  is  a  large  deviation  of  his  ex¬ 
trapolated  results  from  the  actual  pressure  drop  through  the  packed  bed. 

4  5 

In  fact,  at  high  Reynolds  numbers,  10  -  10  ,  it  has  been  shown 

that  the  coefficient  of  drag  is  no  longer  a  sole  function  of  Reynolds 
number,  but  a  function  of  bead  diameter  and  velocity,  other  than  the  ve¬ 
locity  effect  in  the  Reynolds  number.  It  has  graphically  been  illustrated, 
(as  shown  in  Figure  16)  that  irregardless  of  the  flow  area  of  the  packed 
bed,  the  predicted  coefficient  of  drag,  Eq.  (15),  accurately  correlates 
the  measured  friction  factor  (F  )  for  whatever  size  particle  the  bed  is 
composed  of. 

Further,  at  high  Reynolds  number,  the  pressure  gradient  within  a 
packed  bed  is  a  function,  solely,  of  the  average  gas  pressure,  (or  den¬ 
sity),  in  the  bed.  This  pressure  gradient  is  accurately  predicted  using 
Eq.  (16). 

It  must  be  emphasized  that  the  formulated  correlation  for  the  coeffi¬ 
cient  of  drag,  Eq. (15) ,  and  the  pressure  gradient,  Eq.  (16),  were  de¬ 
veloped  for  high  velocity  (and  high  Reynolds  numbers) .  These  correla- 

3 

tions  are  not  applicable  for  a  Reynolds  number  realem  less  than  10  . 
Additional  work  is  now  in  order  to  develop  an  additional  correlation  that 
will  provide  the  transition  equation  from  the  relatively  low  velocity  (and 
lower  Reynolds  number)  conditions,  as  developed  by  either  Ergun  or  Kuo 
and  Nydegger,  to  the  conditions  tested  here. 
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APPENDIX  A 
REDUCTION  FORMULA 


Robbins  and  Gough  [Reference  6]  were  able  to  obtain  pressure  drop 
data,  for  high  pressure  flow  through  packed  beds  of  various  sized  par¬ 
ticles.  Knowing  the  mass  flow  rate,  porosity,  Reynolds  number  and  size 
particle,  they  were  able  to  compute  a  friction  factor  from  the  following 
equation 


F 

s 


(P.2  -  P.2)  <i>3  A2Db  gc 
2(l-<j>)  rfi2  RTL 


(A-l) 


where 


F 

s 


P. 

1 


■pj 


A 


D 


b 


friction  factor 


pressure  at  station  i,  j,  i  >  j 
porosity 


area 

6Vb 

effective  diameter  of  particle  (-? — ) 

6b 


m  =  mass  flow  rate 

R  =  gas  constant 

T  =  temperature 

L  =  length  between  gages  i  and  j 

g  =  constant  to  reconcile  units 

°c 


Multiplying  both  sides  of  (A-l)  by  Re 


where  Re  = 

4> 


P  <j>U  D. 
U 


-  n  <■,  (P.  -  P.)  (P.  +  P.)  <t>  A2  D.  g„ 

.  ReO-il  _  J  i  J  i  P  Sc 

s  .2  L  2  RT  .2 

i>  m 


p  4>U  D, 

m  m  b 


yields, 

(A-2) 


Now,  mass  flow  m  =  p  cf>U  A 
m  m 

Pressure  drops  AP  =  P.  -  P. 

J  i 


(A— 3) 
(A-4) 


APPENDIX  B 


NOZZLE 

The  American  Meter  Sonic  Flow  nozzle  provides  sonic  flow  at  the 

throat  of  the  nozzle,  and  supersonic  flow  at  the  exit  of  the  nozzle,  and 

in  the  test  section  when  empty.  This  fact  was  confirmed  by  running  a  test 

at  400  psig  through  an  empty  test  section  and  recording  the  dynamic  pressures. 

The  nozzle  as  purchased  has  an  exit  diameter  to  throat  diameter  ratio 

Ac 

of  0.6875  to  0.3750  inches,  or  an  area  ratio,  —  =  3.361.  For  ideal  isen- 

A* 

j 

tropic  choked  flow,  Ae/A*  and  Pe/Po  are  related  by 


Ae/A*  = 


Pe 

Po 


Y+l 


where  Y  =  vy 


f 2 

2  (Y-l) 

(Y+lJ 

2  Y 
Y-l 


-  [itr t 


(B-l) 


Solving  Eq.  (B-l)  implicitly,  one  determines  that  Pe/Po  =  0.03916, 
and  hence  Pe  =  16.23  psia. 

We  in  fact  measured  17.62  psia,  16.37  psia,  and  16.37  psia,  at 
gauges  1,  2,  and  3  respectively.  The  differences  in  the  actual  and  ideal 
pressure,  16.23  psia,  were  due  to  the  real  gas  effects  not  being  taken  into 
account.  In  addition,  the  low  gauge  pressures  recorded  on  the  gauges  were 
abou^  2.0  psig,  and  noting  at  ^0.5  psig  accuracy  of  the  gauges,  this  would 
also  account  for  the  small  discrepancies. 

Now  1^-  and  Mach  number  are  related  for  isentropic  flow  by 


61 


Solving  Eq.  (B-2)  for  M,  gives  one  the  Mach  number  in  the  test  section, 

M  =  2.76. 

When  the  test  section  is  packed,  the  flow  can  no  longer  remain  super¬ 
sonic,  but  in  fact  is  reduced  to  subsonic  velocities.  At  a  test  pressure  of 
400  psig,  with  a  packed  test  section  the  average  velocity  in  the  test  section 
is  determined  by  mass  flow,  i.e.. 


.  rc  .  =  U 

o  4>  A  avg 


where 


Pl+P2+P3 


3RT, 


T.S. 


(B-3) 


(B-4) 


giving  an  average  gas  velocity  Uavg  of  124.876  ft/sec. 
and 

U 

M  =  -m — 


/yRT 


CB-5) 


T.S. 


yields  NL,  *  c  ...  =  0.11 

1  Test  Section 

Hence,  the  sonic  nozzle  provides  supersonic  flow  prior  to  the  packed  test 
section,  but  once  the  flow  enters  the  test  section  the  drag  reduces  the 
velocity  to  subsonic  flow. 

The  purpose  of  utilizing  the  supersonic  nozzle  is  to  assure  choked 
flow  through  an  isentropic  contriction  so  that  the  mass  flow  can  easily 
be  calculated  from  the  relation 

m  =  C*  Po  A*  /  /  RT  (B-6) 


where  C*  =  F  •  CD>  and  CD  is  an  effective  discharge  coefficient. 


Equation  (B-6)  includes  a  non-ideal  correction  coefficient,  C*,  to  the 
classical  choked  flow  equation  as  given  in  Reference  10. 


When  calculating  the  mass  flow  through  the  sonic  flow  nozzle,  real 
gas  effects  must  be  taken  into  account.  In  accordance  with  Johnson  [10], 
it  has  been  found  that  real  gas  effects  cause  definite  corrections  in  the 
ideal  gas,  isentropic  flow,  calculations  for  mass  flow  through  a  nozzle. 
Johnson  has  tabulated  a  Critical  Flow  Factor,  C*,  which  takes  into  account 
the  compressibility  and  real  gas  effects  for  mass  flow  through  a  nozzle, 
with  throat  area  A*,  upstream  temperature,  Tq,  and  pressure,  P  ,  and 


R 


R_ 

MW 


53.32 


lbf-ft 

lbm-°R 


for  air. 


A*C*P 

m  =  - 2  (B-7 

/■rt- 

o 


For  ideal  isentropic  flow,  theoretical 

Y+L 


c*  =  r  =  /"y  -At 

'  Y+l 


2(Y-1) 


=  0.6847  if  y  =  1.4 


(B-8 


Johnson  showed  that  the  critical  flow  factor  varies  with  the  plenum 
stagnation  pressure  and  temperature  from  0.684  to  0.770  which  is  a  12% 
variation  from  the  theoretical  ideal  gas  critical  flow  factor  of  0.6847. 


APPENDIX  C 


CALIBRATION 

The  Bourdon  gauges  were  calibrated  by  using  a  Refinery  Supply  Co., 

Model  3S260,  dead  weight  tester.  A  weight  was  applied  to  the  dead 
weight  tester  which  subsequently  applied  a  pressure  to  the  Bourdon  gauge. 

The  respective  gauge  reading  was  recorded.  A  plot  of  the  gauge  reading 
versus  the  applied  pressure  was  prepared  for  each  Bourdon  gauge.  This 
gave  a  calibration  curve  for  each  Bourdon  gauge. 

A  similar  process  was  followed  with  the  Setra  System  Transducers,  how¬ 
ever,  the  output  from  the  gauge  was  recorded  on  a  voltmeter  and  also  an 
oscilloscope  .  Figure  Cl  shows  the  oscilloscope  calibration  display.  Two 
calibration  curves  were  constructed  for  each  transducer,  one  for  the  oscillo 
scope,  and  the  other  for  the  voltmeter.  On  these  curves  the  pressure  was 
plotted  versus  the  voltage. 


Figure  Cl.  Oscilloscope  Calibration  Display 
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APPENDIX  D 
TEST  PROCEDURES 


It  is  imperative  that  the  same  test  procedures  are  followed  for 
every  test,  both  steady  state  and  transient,  so  that  the  resulting  data 
can  be  accurately  compared.  For  the  steady  state  tests,  the  procedure 
is  numerated  below. 

Steady  State 

1.  Close  Grove  Regulator 

2.  Turn  on  main  air  supply 

3.  Turn  on  oscilloscope  (voltmeters)  and  storage  facility 

4.  Power  to  thermocouple  amplifier 

5.  Power  to  pressure  transducer  amplifiers 

6.  Remove  test  section  from  apparatus 

7.  Pack  test  section  with  appropriate  size  beads 

8.  Replace  test  section 

9.  Compress  bed  with  air  pressure  that  is  higher  than  any  test 
pressure  that  will  be  used. 

10.  Adjust  scope  to  zero  output 

11.  Erase  oscilloscope 

12.  Have  oscilloscope  camera  ready 

13.  Open  grove  regulator  to  the  desired  test  pressure 

14.  Run  for  60  seconds.  Recording  pressure  on  Bourdon  gage  2 
and  thermocouples  1,  2  and  3. 

15.  Close  Grove  Regulator 

16.  Inspect  oscilloscope  stored  display  for  any  fluctuation  in 
trace,  indicating  unsteady  flow. 

17.  Take  picture  of  stored  traces  on  oscilloscope 

18.  Erase  and  zero  scope 

19.  Make  run  at  next  higher  pressure,  starting  at  step  11. 

In  the  transient  test,  the  test  apparatus  is  modified  as  explained 
in  Chapter  2.  The  test  procedures  are  listed  below. 


65 


Transient 

1.  Close  Grove  Regulator 

2.  Turn  on  main  air  supply 

3.  Turn  on  oscilloscope  and  activate  storage  facility 

4.  Power  to  thermocouple  amplifier 

5.  Power  to  pressure  transducer  amplifiers 

6.  Remove  test  section  from  apparatus 

7.  Pack  test  section  with  beads 

8.  Replace  test  section 

9.  Compress  bed  with  air  pressure  from  the  9  tank  bank,  same  as 
in  steady  state  case. 

10.  Adjust  scope  to  zero  output 

11.  Erase  oscilloscope  screen 

12.  Have  oscilloscope  camera  ready 

13.  Ensure  ball  valve  is  closed 

14.  Open  compressed  air  bottle  valve 

15.  Activate  sweep  of  oscilloscope  and  open  ball  valve  simultaneously 

16.  Let  air  bottle  empty  until  pressure  in  plenum  drops  below  100 
psig,  then  close  ball  valve. 

17.  Continuously  monitor  temperature  in  plenum  tank. 

18.  Stop  sweep  of  scope 

19.  Take  picture  of  stored  traces  on  oscilloscope 

20.  Erase  and  zero  scope 

21.  Repeat  procedure  with  new  air  bottle  and  a  different  size  bead 
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APPENDIX  E 


TRANSIENT  CHAMBER  EMPTYING 


In  order  to  calculate  the  transient  decay  of  pressure  within  a  tank 

of  gas,  one  must  start  with  a  tank  of  volume  V  ,  initial  pressure  P  ,  anc 

initial  temperature  T  .  The  gas  is  specified  by  the  gas  constant  R  and 

specific  heat  ratio  y.  Assume  that  the  pressure  in  the  tank  is  uniform 

Pa  2  -X. 

(no  waves),  P  <  P  ,  and  <  (— p)  Y-l  (flow  is  choked).  For  choked 

a  o  y  Y+i 


flow. 


P  A*  r 


where  r  =  /~y  (~)  2(y_1)  =  0.6847  (for  air) 


Now,  also 


dt  dt 


(pv)  =  -  v  dp/dt 


Assume  isentropic  expression: 


and  from  (E-2a) 


1-Y 


dM 

dt 


Y  L 


17y 


(p  ) 

v  CT 


]  [p  "  ]  at 


Substituting  equation  (E-l)  for  the  left  hand  side  of 

ill 


( 


a*  r  p 

/  RT 


fv  _ 1 

.  Y  1/Y  . 


Y  d£ 
dt 


Multiplying  (E-6)  by  Pq/Po  and  substituting  equation  (I 
into  (E-6)  yields. 


P  (r//T)  A*  p  [i-(J^)] 

- -  (£-) 

rr~  o 


1-Y 

7  po^  Y 

o 


Substituting  the  state  equation,  PQ  =  P0RT0»  into  CE- 7 

r  a-  /TPT  ■  -  i  (f-)  2Y 

o 


Now  define  P  =  P/PQ  and  rearrange  (E-8) 


dP 

dt 


v/A* 


♦ ; 2Y  =o 


L  y  /W  r  -i 

1  o 


Define  t 


v/A* 


res 


y  r  /RT 

1  n 


(units  of  sec.) 


and 


Thus 


Is  t/Tres 


(nondimensional) 


dP 

dt 


♦  (P) 


3Y-1  q 


wi 


ith  initial  conditions  t  =  0,  P  =  1 


-5)  yields. 


-3),  T/Tq  =  (I 


d(P/PQ) 

__ 


,  one  obtains 
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Solving  equation  (E-10)  yields 


1-n 


-  t  = 


1-n 


(E- 11) 


5y_  i 

where  n  =  -4 —  and  C  =  constant 

2y 

Applying  the  initial  conditions  and  solving  for  the  constant  in  equation 
(E-ll)  gives 


C  = 


CE- 11a) 


Hence 


1 

^  1-n 

P  =  [1  -  (1-n)  t] 


For  y  =  1.4  (air)  equation  (E-12)  becomes 

P  =  [1  ♦  (y  t)]’7 
or  t  -  7  [P  ^  -  IJ 


(E-12) 


(E-13) 
(H -14) 


For  our  specific  case  the  air  bottles  that  are  being  used  have  a 

volume  v  =  309  ff\  P  =  2414.7  psia,  T  =  530°R,  R  =  S3. 3  ,  and 

o  r  o  ’  lbm  R 

Y  =  1.4.  In  order  to  empty  the  gas  bottle  until  conditions  become  un¬ 
choked  (P  =  27.84  psia),  then  P  =  =  .0115,  and  t  =  6.248.  There¬ 

fore,  it  takes  2.73  minutes  (164  seconds).  Figure  E-l  is  a  graph  of 
nondimensional  pressure  versus  nondimensional  time. 


1 

i 

I 


APPENDIX  F 

REDUCTION  OF  EXPERIMENTAL  DATA 


The  experimental  data  was  collected  as  stated  in  Chapter  2.  The 
pressure  transducer  readings  from  the  voltmeters  were  recorded  in  units 
of  volts.  The  transducer  calibration  charts  were  then  used  to  convert 
the  voltmeter  readings  to  pressure.  The  steady  state  pressure  readings 
on  Bourdon  Gauges  1  and  2  were  recorded  on  the  data  sheet,  as  well  as 
the  thermocouple  readings  at  position  1,  2,  and  3.  A  typical  data  sheet 
is  shown  in  Figure  F-l. 


Bourdon 

Reading 

Actual 

Gauge 

PSI 

PSI 

1 

370 

368 

2 

100 

100 

Thermo- 

Couple 


1 

88 

548 

2 

88 

548 

3 

85 

545 

Pressure 

Transduc 


Oscilloscope 

Displacement 


Volts/Div 


m  Volts 


Porosity:  0.3868 

Mass  Flow  Rate:  0.2881  lbm/sec 


Pressure 

Psi 


79.25 

61.5 

43 


Size  Beads:  3mm 
Reynolds  #:  10817.46 


Figure  F-l.  Typical  Two  Phase  Flow  Data  Sheet. 
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After  the  series  of  test  with  the  specific  size  beads,  in  this  case 
3mm,  is  completed,  the  beads  in  the  test  section  are  weighted  to  deter¬ 
mine  their  total  mass.  This  mass  is  then  divided  by  the  density  of  the 
3mm  bead  to  compute  the  porosity  in  the  test  section,  as  shown  below. 
Mass  of  3mm  beads:  657.447  grams. 


Mass  _  657.447  grams 

density  42.292  gram/in* 


15.545  in 


3 


4>  =  1  -  Vol  beads _ 

Vol  Test  Section 


1 


.  .3866 

25.352  in'* 


=  0.3868 


The  mass  flow  is  computed  using  equation  (B-7).  Note  all  calcula¬ 
tions  are  performed  using  absolute  pressure. 

A*C*P 


m  = 


/  RT 


(irin2) (ft2/144  in2) (0.68684) (1 14 . 381bf/in2) (144in2/ft2) (32. 21bm-ft) 
_ _ _ _ _  lbf-sec^ 

/  (53.3  lbf-ft)  (548  °R)  (32.2  lbm-ft  ) 

/  lbm  "®R  ...  2 

lbf-sec 

m  =  0. 2881  lbm/sec 

The  average  density  is  computed  from  equation  (B-4)  and  average  ve¬ 
locity  from  equation  (B-3). 


tPl  1  .  1(79.25  «.<iLS«43)  .  ll38  lbf/i„2J(1441ai1 

3  R  T_  ~  ^ 


T.S. 


ft 


it.  f  /*  ^  r  r-  »  (S 
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p  =  0.3729  lbm/ft'3 

ra 


j  _  _ _  (,2881  lbm/sec) _ __ 

avS  ~  pm  *  Apipe  (.3729  lbm/ft3) (. 3868)  (it  in2)(ft2/144  in2) 


U  =  91.55  ft/sec 
avg 


The  formula  for  Reynolds  number  is  given  by  equation  (4a) 


U  0,  U  \ 

Re  _  o  JL_b.  =  D  -2X& _ ^  =  - 

Re  "  p  u  p  u  A 


D. 

m _  b 

U 


pipe 


Re  = 


(.2881  lbm/sec)  (.01003  ft) 


(tt  in")  (ft“/144  in")  (1.2248  x  10  lbm  ) 

sec-ft 


Re  =  10817.46 


The  reduced  experimental  data  for  a  typical  test  using  a  3mm  and 


0.96mm  bead  in  a  2  inch  test  section  are  tabulated  in  Table  F-l. 


Table  F-l.  Reduced  Experimental  Data  (2  in.  I.D.) 

3mm  Bead  Data,  $  *  0.3868 


AP  1-3  (Psi) 

0  (lbm/ftS) 
m 

A  (±^) 
vsec' 

U  (Fps) 

avg  r 

Re  * 

11 . 25 

0. 1427 

0.1008 

83.71 

3784.93 

19.  25 

0.2201 

0.1649 

88.78 

6191.81 

36.25 

0.3729 

0.2881 

91.55 

10817.46 

32.  2S 

0.53S4 

0.4138 

91.58 

15536.61 

67.75 

0.6966 

0.5409 

92.01 

86.50 

0.8593 

0.6685 

92.18 

24851.42 

102.75 

1.0235 

0.7944 

91.97 

29532.02 

111.0 

1 . 1967 

0.9311 

92.20 

34614.14 

131.5 

1.3451 

1.0592 

93.31 

38991 . 96 

0.96mm  Bead  Data,  $  *  0.3948 


Plenum 

Pressure  AP  1-3  (Psi) 


0m  Clbm/ft3)  m  O 


0.1401 


0.0267 


0.3979 


0.5803 


0.7497 


0.9206 


1.0915 


6560.98 


8092.35 


9558.90 


10939.84 


400 


1.4294 


1.0737 


87.22 


12445.36 
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